Lorentz-violating effects on topological defects generated by two real scalar fields 
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The influence of a Lorentz-violation on soliton solutions generated by a system of two coupled 
scalar fields is investigated. Lorentz violation is induced by a fixed tensor coefficient that couples 
the two fields. The BogomoPnyi method is applied and first-order differential equations are obtained 
whose solutions minimize energy and are also solutions of the equations of motion. The analysis 
of the solutions in phase space shows how the stability is modified with the Lorentz violation. 
It is shown explicitly that the solutions preserve linear stability despite the presence of Lorentz 
violation. Considering Lorentz violation as a small perturbation, an analytical method is employed 
to yield analytical solutions. 
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I. INTRODUCTION 



Since the demonstration that spontaneous breaking of Lorentz symmetry is allowed in the high energy context of 
string theories [l[ , Lorentz-violating theories have been extensively studied in diverse low energy systems and used as 
an effective probe to test the limits of Lorentz covariance with direct consequences on the Planck scale physics. The 
great majority of such investigations take place in the framework of the Extended Standard Model (ESM), conceived 
by Colladay and Kosteletcky [2[ as a extension of the minimal Standard Model of the fundamental interactions. The 
ESM admits Lorentz and CPT violation in all sectors of interactions by incorporating tensor terms (generated possibly 
as vacuum expectation values of a more fundamental theory) that account for such a breaking. Actually, the ESM 
model sets out as an effective model that keeps unaffected the SU(3) x SU(2) x U(l) gauge structure of the underlying 
fundamental theory while it breaks Lorentz symmetry at the particle frame. There is a large amount of work in the 
literature [H, 0| searching for experimental and observational evidence of Lorentz breaking, in an attempt to put lower 
bound on the parameters. 

Topological defects are actually one of the most active research branches of theoretical physics, with profound 
connections with condensed matter physics, gravitation and field theory [5]. In particular, kink solutions are the 
simplest topological defects which appear after symmetry breaking of models with one or more scalar fields. Explicit 
solutions have been worked in several models involving one or more scalar fields @, including gravity Q, nested 
defects Q, junctions H or connections with other areas of physics [Io| . 

The influence of Lorentz violation on kink solutions was investigated in ref. [ll| , where the first-order framework 
was applied for some class of models finding explicit solutions in some interesting cases. Here we pursue this problem 
studying another class of models with some new results including explicit solutions, phase space analysis, stability 
analysis and perturbative solutions. This paper is outlined as follows. In Sec. II, we present our model with Lorentz 
violation and we apply the Bogomol'nyi method to find first-order equations that minimize the energy and solve the 
equations of motion. The behavior of numerical solutions are studied in the phase space for several values of the 
Lorentz breaking parameter. In Sec. Ill we prove that our solutions are stable and we find that the model support 
traveling waves. Further, in Sec. IV we find explicit perturbative solutions for the fields. Our final conclusions are 
presented in Sec. V. 



II. THE MODEL AND FIRST-ORDER EQUATIONS 



We consider a class of models with two real scalar fields 



C = \d^d^ + \d» X dvX + k^d^dpx - V(tf>, x), (1) 
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where the coupling tensor characterizes Lorentz violation. The tensor coefficient is a 2 x 2 matrix, 



f3 a 
a (3 



and is an element that couples the scalar fields, introducing Lorentz violation and inducing an asymmetry in the 
propagation of the traveling waves present in these models. The Euler-Lagrange equations can be written as 

t - 4>" + P(X + x") + 2«X + ^ = 0, (2a) 
x -x" +l 3($ + <j>") + 2a<fi' + V x = 0, (2b) 

where = dV/dcf), V x = dV/dx, and prime and dot stand for space and time derivative, respectively. In the case of 
static solutions, we get 

-0" + p X " + V 4 = 0, (3a) 
-x" + H" + v x = 0. (3b) 

The energy- momentum tensor, 

Tiiv = 1 [g^vf + d» X d v X ~ W a d aX d v <j> + k a ^d a (j)d v x ~ 9^C] , (4) 
provides the energy density of the static solutions: 

t 00 = \ + x' 2 - m'x! + W{4>, X)] ■ (5) 

A useful form for representing the potential V((f>, X ) can be attained from Eqs. (|3a|) and (|3b[) . as we get 

4z(-y 2 -lx' 2 + Wx' + v) =0. (6) 



dx \ 2 2 
Then a potential satisfying 

V = + \x' 2 ~ P<f>'x' (7) 

would lead to the second-order equations of motion for the fields </> and x- This shows that the gradient and potential 
portions of the energy do not contribute equally, the difference being proportional to the Lorentz violating coefficient 
/?• 

In order to obtain first-order differential equations whose solutions are also solutions of the equations of motion we 
use the Bogomol'nyi approach [12j, looking at the energy density T 00 and trying to extract simple conditions that 
minimize the total energy. We introduce a smooth function W = W{4>, x) an d consider the ansatz 

4>' = W^ + sxx 1 , (8a) 
X' = W X + s 2 <j>'. (8b) 

Substituting this into Eqs. ([5]) and ([7]) leads to 



7 = lw 2 + x' 2 )-Wx' + lm + Sl x') 2 + l(w x + sri') 2 >'"Y 



= + x'W X ) + <f>'x'(si + S2- 2/3) + \ (</>'- - s lX ') 2 + \{x'~W x - S2 cp') 2 . 

For Si + S2 = 2/3 and when the fields <f> and x satisfy the ansatz given by Eq. ©, we have T 00 = e = dW/dx, and the 
total energy E is minimized to the value 

/>oo 

E= / e(x)dx= \AW\. (9) 
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In the remaining of this section we take si = s 2 = p. Explicitly, the ansatz given by Eq. (|5J) results in two first-order 
differential equations that minimize the energy 



X —(W^ + f3W x ) (10a) 



Note that Eq. ([7]) can be written as 



1-p 

x = y^{w x + PW^) (iob) 



I , I 

_ 

1 -/3 2 V2 



V{4>,x) = r r ^[^ + 7 ffi+PW*W x ) i (II) 



and for f3 2 < 1, which concur with the fact that any Lorentz violating parameter should be very small, the potential 
is positive. The simplest solutions for the equation (|10[) are the homogeneous solution. These solution can be found 
using the W$ — W x = that are minima of potential. 



A. An Example 



Now we will consider a particular potential. We choose the function W introduced in Ref. 
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corresponding to the potential 

V{cj>, X )-- 



1 ~f3 2 



W r W,x) = 0-§0 3 -r^c a 



(1-4,'- rx'Y + 2rW - 2r/3fo(l - 4>~ rtf) 



(12) 



(13) 



This potential depends explicitly on (3, but the minima are (4>,x) = (0, i-v/l/r) and (4>iX) — (^IjO) an d they do 
not depend of (3. These points correspond to the following values of W: Pi = (4>, X ) = (— 1, 0) =>■ W\ = —2/3; 
P 2 = (1,0) =► W 2 = 2/3; P 3 = (0, xfTfr) W 3 = 0; P 4 = (0,-y/TJr) =4> W 4 = 0. In this way we can possibly have 
five BPS sectors connecting Pi «-» P 2 , -Pi *H> -P3, -Pi <-> -Pi, P2 <-> -Pi, -P2 <-> -Pi and one non-BPS sector connecting 
P 3 f+ P 4 . 

Now we analyze Eqs. (|10l) on the (0, x) plane, looking for explicit solutions. Those equations reduce, in the model 
considered here, to the following first-order equations for the fields: 



1 



1 -/3 2 
1 



[1 - 4> 2 - rx 2 - 2/3r<p X ] 
[-2r<f>x + ^l-^-rx 2 )} 



The case (3 = lead to 



<t>' = l-tf-rx 2 , 
X 1 = -2r0x, 



which can be integrated to give the general orbit constraint 

,0 r 



2r-l 



X 1 + C X ~ + 1 



(14a) 
(14b) 



(15a) 
(15b) 



(16) 



A special case is the elliptic orbit for C = 0, with <fr 2 + x ^ 2r X 2 — !• As this particular orbit connects the minimum 
energy points Pi and P 2 , it is of particular interest with explicit solutions 



(x) = tanh(2rir), x( x ) = \l 2 scch(2rx) 



(17) 



Another special orbit is the line x = (C — > 00). In this case, we take <j>(x) = tanh(x) 
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FIG. 1: Profile of the orbits in the (<j>, \) plane, given by Eqs. (|14|l . for ft = 1/3 and r = 1/3. The orbit in black (dot-dashed 
line) corresponds to odd solutions, which have the symmetry <j>{— x) = —(j>(x) and x{~ x ) = ~x( x )- 



We have been unable to find explicit solutions for general nonvanishing values of p. However, numerically we could 
obtain orbits and the corresponding solutions (f>(x) and x( x )- I n Fig.[U we plot the orbit profile of the Eqs. (fT4"|) . The 
green (solid) lines connect the minima Pi and P%. The blue (dashed) lines connect the minima Pi and P2 to P3 and 
P4. Note that, in contrast with the case of vanishing /?, there is no solution with one constant field. 

We plot three pair of solutions in the next three figures, and in Fig. [2] we show that the <j>(x) (left panel) and x( x ) 
(center panel) solutions are related to the central orbit highlighted in black in Fig. [TJ which connects PI and P2. 
The right panel shows the energy density, and we note that there are two small energy peaks in addition to the large 
central peak. These solutions have energy Eb =4/3. In Fig. [31 we show the solutions and the respective energy 
density related to the orbit connecting PI and P3. The <j> solution is not a monotonic function anymore, with a 
minimum at x ~ 3.5. Finally, in Fig. [U we show the solutions that connect P2 and P3. We note that, in contrast 
with the P = case, the BPS solutions that connects PI to P3 do not have the same form as the BPS solution that 
connects P2 to P3, as we can see in Fig.[3]and in Fig. [4] Nevertheless, both solutions have the same energy Eb — 2/3. 



FIG. 2: The <f>(x) (left panel) and x( x ) (right panel) solutions and energy density (right panel), for ,5 = 1/3 and r = 1/3. These 
solutions connect the points PI and P2. 



In the case of /3 being very small, we can find approximate results. However , we will present a perturbative analysis 
for $ < 1 in Sect. [TV] 



In order to analyze the stability of the solutions, one should start with the full partial differential equations (|2aj) 
and (l2b| . For a — 0, for simplicity, we get to 



It is interesting to observe that the Lorentz violation does not jeopardize the linear stability of the differential equations, 
which turns stable the associated kinks solutions. Such demonstration is performed according to the usual way. With 




III. STABILITY ANALYSIS 



^ - 4>" + P(X + X") + dV/dcP = 0, 
x - x" + P(4> + 0") + dV/dx = 0, 



(18a) 
(18b) 



this purpose, we consider general fluctuations, r/, £, on the field configurations: 

<p(x, t) = (j>(x) + r)(x, t), x(x, t) = x(x) + £(x, t). 
With this, the differential equations (fT8l) take the form 



a 2 ij - ot\r)" + V Vlipi J] + V VlV>2 £ = 



where a>i = (1 — /3), a 2 = (1 + 0) and 



(19) 

(20a) 
(20b) 

(21) 
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FIG. 3: The same as in Fig. [2] for solutions that connect the points PI and P3. 




FIG. 4: The same as in Fig. [2] for solutions that connect P2 and P3. 
Such system may be written in a matrix form: we use r)(x, t) = rj(x) cos(wi) and t) = £(x) cos(oji) to get 



H 



a 2 ij 



where 



H 



d 

vin-* inr, 



'<Pl<t 



-OL2 



Finally, it is possible to show that H = S^S, with: 

S 



PiW Vl , 



P2W V1V2 



(22) 



(23) 



(24) 



with pi — a-Jax, p 2 = b^/a^. Eq. (|22p is not an eigenvalue equation but its structure is enough to demonstrate 
stability. Indeed, such equation may be written as H\*f?) = u;^M|\f r ), where M is a diagonal 2x2 matrix with 
Mn = a 2l M 22 = ai. Considering the internal product = <^ r |S' t S , |^ r ) = IS*!*)! 2 , we obtain IS"!*)! 2 = 
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w n( a 2M 2 + Q!i ICI 2 ) j which assures positivity of w 2 , since a\,a 2 > 0. The positivity of w 2 implies real values for u, 
which yields well behaved solutions for any time, meaning stability. 

Another study concerns the propagation of traveling waves in this model of two coupled fields. First of all, we note 
from Eqs. (JTSJl that 

a2(f>i — aiip" + dV/dipi = 0, (25a) 
aiv3 2 - a 2 <f2 + dV/d<p 2 = 0, (25b) 

If such a model entails the presence of traveling waves, then there should exist moving waves endowed with the same 
form as the static solutions, which satisfy a\Lp'{ = dV/dtpx, ot2ip 2 ' = dV/dip 2 - In order to properly examine this issue, 
the fields are written in the form tpi = ipx(ux),<p2 = ^2(^2), with u% = Ji(x — v\t),U2 — 72(2^ ~ ^it), where 71,72 
being constants and v\, and v 2 represent the propagation velocities. Replacing ^i(wi), ^2(^2) into eqs. (|25[) . one gets 

aid 2 (pi/du 2 = dV/dcpi, (26a) 

a 2 d 2 ip 2 /du 2 = dV/dip 2l (26b) 

after making the following identification 

71 - [1 - (a2/a 1 )v 2 ]- 1 ' 2 , 72 = [1 - (a 1 /a 2 )v 2 ]- 1 ^. (27) 

The attainment of Eqs. (|26al26bl) confirms that this model supports traveling waves. The phase velocity (v p h — ui/k) 
of such waves is given by the dispersion relations coming from the free version of Eqs. (125al25b[) . —012k 2 + ol\uj 2 = 0, 
— a.\k 2 + «2W 2 = 0, which implies v\ = a 2 / a\, V2 = ±y ai/«2- Considering the redefinitions ([21]) . the stationary 
wave configurations associated with the fields <j), \ are composed by two travelling waves of velocities v\ , V2 ■ It is 
interesting to note that i; 2 < 1 while vi > 1, and that this outcome prevails for the group velocity as well, with 
v g2 < 1 while v gi > 1. 



IV. ANALYTICAL METHOD FOR PERTURBED SOLUTIONS 



Taking into account that Lorentz violation is a small effect, a perturbative approach is well suited to deal with 
it. In this sense, an analytical method may be used to construct explicit first-order solutions for the differential 
equations (|18[) constrained to the situation where the LV parameter /3 is small (since its presence is addressed as a 
perturbation on the Lorentz symmetric case). In order to obtain first order analytical solutions, we should first define 
the unperturbed system (/? = 0) as 

4> -4>" + dV/dcf) = 0, (28a) 
X - X" + dV/dx - 0, (28b) 

which provides unperturbed solutions (labeled as $0 and xo), whose form depend on the adopted potential V. In this 
section we consider V(<j), x) — Vi{4>) + ^(x), which makes the above equations uncoupled. 
In the static perturbed case, the differential equations are: 

-0" + fix" + dV/d<j) = 0, (29a) 
- X " + H" + dV/d X = 0. (29b) 

Note that these equations are coupled only by /3-dependent terms. In order to solve it, we propose the following 
ansatz 

<f>{x) =Mx) + Pg(x), x(x)=Xo(x) + Pf(x), (30) 
where g and / are functions to be determined. Replacing such ansatz in eqs. (|29al29b|) , it results 

-9" + Xv+9d 2 V/d(j>% =0, (31a) 
-f" + ^ + fd 2 V/d X 2 \ X0 =0, (31b) 

where 0o and xo are the corresponding unperturbed solutions. These are second-order differential equations, whose 
solutions may be achieved in an analytical form for certain potentials. Aa a first example, let us consider the potential 
such that 

Vi(0) = (l-0 2 ) 2 /2, F2(x)-(1-X 2 ) 2 A (32a) 
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which provides the well-known solutions, 4>q(x) = ±tanh(x), Xo( x ) = ±tanh(x), besides the homogeneous ones: 
<t>o{x) = Xo( x ) — 0- I n this way, Eqs. (|31a[) - (|31bl) give 

-9" + Xo + ff(-2 + 6^) = 0, (33a) 
-/" + 0o+/(-2 + 6Xo)=O. (33b) 

A first situation consists in choosing <j>o(x) — ±tanh(x), Xo( x ) = il> which leads to the system of equations 

-p" + ,g(6 tanh 2 a; - 2) = 0, (34a) 
-/"=F2tanhxsech 2 a; + 4/ = 0. (34b) 

Under the boundary conditions (dg/dx)\ x ^ f00 — (df /dx)\ x ^ f00 — 0, g(0) — a, /(0) = 6, such system supports an exact 
analytical solution given as follows 

g(x) = a sech 2 x, (35a) 
f{x) = be" 2x ± 2[1 - ln(2coshx)] sinh2x ± 2x cosh 2x =p tanh x(l + 2cosh2x). (35b) 

With this, the classical solutions for the Lorentz-violating system are determined at first order in accordance with 
eqs. (j30]l. 

A second case is </>o(x) = ±tanh(x), Xo( a; ) = tanh(x) with the conditions {dg / dx)^^,^ = (d//dx)| x _ > . 00 = 0, 
g(Q) = a, /(0) = b, which gives 

g(x) = (a + x/2) sech 2 i, /(x) = (b ± x/2) sech 2 x. (36) 

Once we have obtained first-order analytical solutions, it is also possible to carry out the corresponding first-order 
energy corrections for the solutions already known. The starting point is the energy density given by eq. ([5]), written 
for eqs. ([3D]), which implies T 00 = T 00 + T^°, where: 

T°° = iK 2 + Xo 2 + 2n0,x)] (37a) 
T°° = ^g 1 + x'of - 0oXo + V<f>9 + V x f}. (37b) 

Here, Tq° is the usual energy density (without Lorentz- violation) and Tp is the first-order correction due to Lorentz 
violation. Now, using the potentials of Eqs. (132ap a- (|32al) b. the energy correction takes the form 

T°° = p\<M + x'of ~ 0oXo " 20o(l - 4>l)9 ~ 2xo(l " Xo)/]- (38) 

A first case whose energy is to be analyzed is the one in which <fio(x) = ±tanh(x), xo( x ) = ±1. In this case, the 
corresponding solutions are given by eqs. (|35al35b[) . which lead to the following result 

T 00 = sech 2 x, (39a) 

pOC 

( /3 



T°° = ^(3[6a tanh x sech 4 x] . (39b) 



The energy associated with this density, E = Ta°dx, is zero, since one has to integrate an odd function. Hence, 

the energy is simply the one of the usual case (without LV): E = T§°dx = 4/3. This means that the energy does 
not change up to first-order in /3. 

The second case is given by 4>a(x) = ±tanh(x) and Xo( x ) = tanh(x), whose solutions are given by eqs. (??|3S|). In 
this case, one gets 

T 00 = 2 sech 2 x, (40a) 
= =f y [x tanh x sech 4 x] . (40b) 

The associated energy correction is then equal to =f4/3/3, with total energy being given by E = 4(2 =p /3)/3. Here the 
energy receives first-order correction in /3. 
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V. ENDING COMMENTS 

In this work, we have studied the effects of a Lorentz- violating tensor term on the topological defects associated with 
models described by two real scalar fields. We applied the Bogomol'nyi approach to successfully attain first-order 
differential equations whose solutions minimize the energy and solve the equations of motion. This simplification 
is extremely useful to study numerically the solutions in phase space. We studied a model where the potential 
describing the two scalar fields presents nonlinear couplings between the fields, with a rich structure of minima. As 
we have shown, for general values of the parameter which implements the Lorentz breaking, a numerical investigation 
was performed, and we have found numerical solutions connecting minima of the potential. For general models, we 
performed a stability analysis to show that the solutions of the equations of motion are well behaved in time. Also, we 
have shown that the models support traveling waves in general. Since we where not able to find explicit solutions for 
general values of /3, we studied the case with f3 <C 1, which concur with the fact that the Lorentz breaking should be 
very small. In this case, we proposed an ansatz for the perturbed solutions, finding explicit solutions for the special 
case where the potential is described by the two fields <f) and X, that do not interact with each other. 
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